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Abstract: The Wireless Personal Area Network (WPAN) is one of the fledging paradigms 
that the next generation of wireless systems is sprouting towards. Among them, a more 
specific category is the Wireless Body Area Network (WBAN) used for health monitoring. 
On the other hand, Ultra- Wideband (UWB) comes with a number of desirable features at 
the physical layer for wireless communications. One big challenge in adoption of UWB in 
WBAN is the fact that signals get attenuated exponentially. Due to the intrinsic structural 
complexity in human body, electromagnetic waves show a profound variation during 
propagation through it. The reflection and transmission coefficients of human body are 
highly dependent upon the dielectric constants as well as upon the frequency. The 
difference in structural materials such as fat, muscles and blood essentially makes 
electromagnetic wave attenuation to be different along the way. Thus, a complete 
characterization of body channel is a challenging task. The connection between attenuation 
and frequency of the signal makes the investigation of UWB in WBAN an interesting 
proposition. In this paper, we study analytically the impact of body channels on 
electromagnetic signal propagation with reference to UWB. In the process, scattering, 
reflectivity and transmitivity have been addressed with analysis of approximate layer- wise 
modeling, and with numerical depictions. Pulses with Gaussian profile have been 
employed in our analysis. It shows that, under reasonable practical approximations, the 
human body channel can be modeled in layers so as to have the effects of total reflections 
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or total transmissions in certain frequency bands. This could help decide such design issues 
as antenna characteristics of implant devices for WBAN employing UWB. 

Keywords: reflection; transmission; scattering; propagator; attenuation 



1. Introduction 

Wireless Body Area Networks (WBANs) have attracted interest in recent years because of a 
number of promising applications — specifically, in the field of health monitoring. Like everyday attire, 
in a WBAN, several small nodes are placed directly in, on or around the human body. Since WBAN 
nodes acquire their power from rechargeable batteries or by energy harvesting, it is essential that they 
be extremely energy-efficient [1]. Above and beyond the energy efficiency, the nodes are meant to be 
of low complexity to keep costs down, among other things. On the other hand, Ultra -Wideband 
(UWB) communication is a transmission technology that comes with such promises as low-power 
consumption [2], interference robustness [3], high local capacity [4], and less complex hardware, most 
of which are highly desirable for WBANs [5]. One key concern in this regard is about signal 
attenuation which occurs exponentially with frequency. This leads to the need to study electromagnetic 
propagation across human body as medium, with consideration for UWB signal as it relates to 
communication system parameters for implant devices. Particularly, Impulse-Radio (IR) [6] 
transmission appears to be well suited to reduce complexity, since major parts of narrowband 
communication systems such as mixers, RF (Radio Frequency) oscillators, or Phase-Locked Loops 
(PLLs) can be omitted in IR systems [7]. In order to accomplish the requirements mentioned above as 
they relate to energy efficiency and complexity reduction, the distinct behavior of the propagation 
channel has to be taken into account. For WBANs, this has to do with identification of the effects of 
propagation on or around the body. Unlike conventional wireless channels, a human body is rather 
complex in structure. Electromagnetic waves show a profound variation during propagation through 
human body as the reflection and transmission coefficients are highly dependent upon the dielectric 
constants [10], in addition to the frequency. Therefore, it is a challenging task to make a complete 
characterization of the human body channel. Due to the differences in structural materials such as fat, 
muscles, blood etc., electromagnetic wave attenuation is different across the different parts of the 
body. The higher the frequency, the more attenuation takes place, which limits the use of high 
frequency or UWB in WBAN. In this manuscript, we analyze the impact of the body channel on the 
signals in different frequency bands. Scattering, reflectivity, and transmitivity have been studied with 
analysis of approximate layer-wise modeling, and with relevant numerical rendering. Pulses, having 
Gaussian profile, have been used in our analysis. We illustrate that the body channel can be 
mathematically modeled as composed of layers, with total reflection and total transmission of the 
signal in certain frequency bands, so as to approximate the propagation effect. The rest of this paper is 
organized as follows: the whole mathematical model for wave propagation in biological media, 
scattering, reflection, and anti-reflection by single layer is given in Section 2; Section 3 describes 
numerical results and concluding remarks are given in Section 4. 
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2. Mathematical Model 



2.1. Wave Propagation through Biological Media 



From the appendix A, if the incident wave is a linearly polarized uniform plane wave travelling 
along the z-direction, then, for E and H, Equation (A. 9) is of the form: 



H = H i e- az e M - pz) i y 

where E i =rjH i 

The intrinsic impedance of biological material n is given by: 
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The Pointing Vector, that is, the power flowing per unit area of cross section (W/m ), gives the 
power density associated with an EM wave: 



P =E : xH ; 



For a uniform plane wave, time-average power flow is given by: 

\E.\ 2 

P = 



— - = -tj\H.\ 



(4) 



(5) 



The permittivity and frequency may also determine how far the EM wave penetrates into the body. 
The term depth of penetration (D p ) usually quantifies this. For objects with homogeneous properties 
and with RFR incident at right angles to the surface, depth of penetration is defined as the distance at 
which the power density is decreased by absorption to about 0.13534 of the body's surface value. 
However, the magnitude of the electric and the magnetic field reduces by a factor of 0.36788. Depth of 
penetration is defined as: 

1 



D„ = 



a 



(6) 



where a is the attenuation constant of the material in nepers per meter. 
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Figure 1. Power absorption in muscle as a function of depth at different frequencies. 
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Figure 2. Variation of Penetration depth with frequency. 
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2.2. Scattering, Reflection, and Anti-Reflection 



Sensor nodes find a human body, when they are placed, to be layered media. Fat, muscles etc. are 
such independent layers. Signals from wireless sensors in human body experience scattering, 
reflection, and diffraction by those layers. In Appendix B, we started with acoustic wave for 
mathematical formulation and then used these formulae for electromagnetic wave. We have 
established the basic equations for the wave propagation through the layered media. Now we will use 



Sensors 2010, 10 



5507 



those mathematical equations for calculating some wave properties such as scattering, reflection, and 
transmission while passing through layered media. 



2.2.1. Scattering by a Single Interface 



In this section we consider the case in which two homogeneous half-spaces are separated by an 
interface at z = 0 (Figure 3): 

IPi l f z > 0 (Kiif z> 0 

The goal of this section is to analyze the scattering problem in terms of right- and left-going modes 
(Appendix B). 

We introduce the local velocities c ; = jKj / pj and impedances ^ = JKj pj and the right- and 
left-going modes defined by: 

\ A Q (t,z) = ^ 1/2 p(t,z) + C 0 /2 u(t,z) 
[B Q (t,z) = -<;- 1/2 p(t,z) + ( 0 1/2 u(t,z) 

\ A 1 {t,z) = (- 1/2 p(t,z)+tf /2 u(t,z) 
[B x {t,z) = -^ 1/2 p(t,z) + £ /2 u(t,z) 

For j = 0, 1, the pairs (Aj, Bj) satisfy the following system in their respective half-spaces: 



z < 0 



z > 0 



(7) 
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which means that Aj (t, z) is a function oft — z/ Cj only, and Bj (t, z) is a function of t + z/c ; only. 

Figure 3. Scattering of a pulse by an interface. 
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We assume that a right-going wave with the time profile / is incoming from the left and is partly 
reflected by the interface. We also assume a radiation condition in the right half-space so that no wave 
is coming from the right. Assume that/is completely supported in (0, oo). We next introduce two ways 
to define proper boundary conditions: 

(I) We can consider an initial value problem with initial conditions given at some time t 0 < 0 by: 

l i xn 
u(t = t 0 , z) = —172 /(t 0 - z/c 0 ), p(t = t 0 , z) = -^—f(t 0 - z/c 0 ) (10) 

As shown in the Appendix B, these initial conditions generate a pure right-going wave whose 
support at t = to is in the interval z G (— oo, c 0 t 0 ), which lies in the left half-space. 
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(II) We can consider a point source located at some point z 0 < 0 and generating a forcing term of 
the from: 



F(t,z) = ( 0 1/2 /(t - z 0 /c 0 )5(z - zq) 



(ID 



As seen in the Appendix B, this point source generates two waves. The left-going wave is 
propagating into the negative z-direction and will never interact with the interface, so we will ignore it. 
The right-going wave first propagates in the homogeneous left half- space and it eventually interacts 
with the interface z = 0. 

In terms of the right- and left-going waves, these two formulations give the same descriptions. We 
have A 0 (t, z) = f(t — z/ c 0 ) for z < 0, and Bi(t, z) = 0 for z > 0, and consequently, at the interface 
z = 0: 



A 0 (t,0) = /(t), 5 1 (t,0) = 0 



(12) 



Note that the delays introduced in the initial conditions (10) and in the forcing term (11) have been 
chosen so that the boundary conditions (12) have a very simple form. 

The pairs (Aq, Bo) and (A\, B\) are coupled by the jump conditions at z = 0 corresponding to the 
continuity of the velocity and pressure fields: 

_ 1/2 (A 0 (t,0)+B 0 (t, 0)\ _ 1/2 (A^t, 0) + g 1 (t,0)\ 



u(t.0) =Q 



p(t,0) = Co 



_ 1/2 ^ 0 (t,0)-g 0 (t,0)^ _ 1/2 /^Qj-giOQ) 



which gives: 



^4i(t,0) 
Bi(t, 0) 



^o(t,0) 
S 0 (t,0) 



^,(+) ) 

2 



(13) 



with r (±) = i (Vd/^o ± VCo/Ci)- Note that - (r (_) ) =1. The matrix /can be interpreted as 

a propagator, since it "propagates" the right- and left-going modes from the left side of the interface to 
the right side. Such a propagator matrix will be called interface propagator in the following. 
Taking into account the boundary conditions (12) yields: 



[ M o 0) H 



lB 0 (t, 0) 



and solving this equation gives: 

B 0 (t, 0) = Xf(t), A 1 (t, 0) = Tf(t) 
where !R and T are the reflection and transmission coefficients of the interface: 



31 = - 



T = 



r« ( 0 + Ci ' rC+) 
These coefficients satisfy the energy-conservation relation: 

Jl 2 +T 2 = 1 
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meaning that the sum of the energies of the reflected and transmitted waves is equal to the energy of 
the incoming waves. Finally, the complete solution for z < 0 in terms of the right- and left- going 
modes is: 

A 0 (t,z) = fit - z/c Q ), S 0 (t,z) = #/(t + z/cq) 

and for z > 0: 

A^t.z) = Tfit-z/ctXB^t.z) = 0 
Using (7-8) we can obtain the pressure and velocity fields (Figure 4). 

Figure 4. Scattering of a pulse by an interface separating two homogeneous half-spaces 
(co, Co, z < 0) and (ci, Cu z > 0). Here the incoming right-going wave has a Gaussian profile, 
Co = Co = h and c\ = Ci = 2. The spatial profiles of the velocity field (a) and of the pressure 
field (b) are plotted at times t = -4, t = -3,..., t = 6. 
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2.2.2. Single-Layer Case: Scattering 

In this section, we consider the case of a homogeneous slab with thickness L embedded between 
two homogeneous half-spaces (Figure 5). Three regions can be described as follows: 

f p 0 if z < 0, ( K 0 ifz<0 

P(z) = \ p 1 if z E [0, 1], K{z) = \K 1 ifzE [0, L] 
{ p 2 if z < 0, { K 2 if z < 0 

We introduce the local velocities c ; = ^K~Jp~ and impedances <" ; = yjKjPj and the local right- and 
left-going modes defined by: 

Aj(t,z) = (- 1/2 p(t,z) + Zj /2 u(t,z), Bj(t,z) = -(- 1/2 p(t,z) + ^ /2 u{t,z) 

with j = 0 for z < 0,y = 1 for z G [0, L], and j = 2 for z = L. The boundary conditions correspond to an 
impinging pulse at the interface z = 0 and a radiation condition at z = L2: 

A 0 (t,0)=f(t),B 2 (t,L) = 0 
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The propagation equations (9) in each homogeneous region show that Aj is a function of t — z/cj 
only and R- is a function of t + z/cj only. The waves inside the slab [0, L] are therefore of the form: 

A-iit, z) = a.i(t — z/c-i), Bx{t,z) = b x (t + z/c{) 
while the reflected wave for z <0 is of the form: 

S 0 (t,z) = b o (t + z/c 0 ) 
and the transmitted wave for z > L is of the form: 

z -Ln 



A 2 (t,z) = a 2 (t 



We want to indentify the functions bo and «2, which give the shapes of the reflected and 
transmitted waves. 



Figure 5. Scattering of a pulse by a single layer. 
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2.2.3. Single-Layer Case: Reflection and Transmission Coefficients 



The unknown functions bo and a2 can be obtained from the continuity conditions for the velocity 
and pressure at the two interfaces. At z = 0, we have: 



A t (t, 0) 




A(t,o)" 


i?i(t,0). 


= /o 


.£ 0 (t,o). 



Jo = 



with r 0 (±) = - (JtJTo ± V?o7^). Similarly, at z = L: 



L5 2 (t,L)J 



^i(t,L) 
LSi(t,L)J 



, Ji = 



.(+) „(-) 
o 'o 

.(-) „(+) 



.(+) „(-) 
l 'i 

.(-) „(+) 



with r^ +>> = i (V^/^i ± V^i/^)- We can write these relations in terms of the functions ay, as: 



"aiCt) 




r/coi 




= /o 


U»o(OJ 



0 



^(t + L/c^J 



which can be solved to get the reflected and transmitted waves. The situation is more complicated than 
in the case of a single interface, because of the time delays ±Llc\. A convenient and general way to 
handle these delays is by going to the frequency domain, so that the time shifts are replaced by phase 
factors. The Fourier transforms of the modes are defined by: 
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a ; (o)) = j aj(t)e ia>t dt, bj{a)) = j aj{t)e ia>t dt 
They satisfy the interface conditions: 



Si (<u) 



= /o 



/(G)) 



a 1 {oS)e 



.Oil 



.d)L 
-I 



(14) 



where we have used the identity: 

j a x {t - L I 'c{)e ib)t dt = j a x {s)e 
Introducing the frequency-dependent matrix: 



/ L \ .a)L 
>(sH — , „ , N i — 

v c i J ds = d^a^e c i 



Jd<o) = 



(+} I (-) -I 

?y e c i J e c i 

, „ .a>L .coL 

(-) ! tt~ (+) 



The second equation of (14) can be written as: 

=JM 



d x (o>) 



(15) 



The syplectic matrix / 1 (a>) is a propagator in the frequency domain. It propagates the right- and 
left- going modes from the right side of the interface 0 to the right side of the interface 1 , and it depends 
on the layer thickness L. Finally, combining the first equation of (14) and (15), we obtain the relation: 



0 

where the frequency-dependent syplectic matrix: 

K 0 (oj) =7i(o))/ 0 = 



fcoO"). 



(16) 



t/(o0 V((o) 



|y(«0 Ufa)! 

is the overall propagator of the slab. Equation (16) shows that K 0 (o>) propagates the right- and 
left-going modes from the left side of the interface 0 to the right side of the interface 1. We 
find explicitly: 



U{co) = r 0 (+) r 1 (+ V c i + r 0 l >e 

.o)L .uiL 



.coL 



(-),.(-)-* 



.<dL 



By solving equation (16), whose unknowns are d 2 {o)) and S 0 (o>) and using the expressions of r. 
we obtain: 

b 0 {do) = a 2 (o)) = f {(d) f {(d) 

where the frequency-dependent reflection and transmission coefficients are: 

.<dL 



(+) 



R{(D) = 



0(co) 



/? 1 e 2 ' c l+/? 0 
2i— 

l+/? 0 /?ie ' c l 



(17) 
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. fflL 

l+/? 0 /?ie c l 

using that |t/(»| 2 - |K(co)| 2 = 1. Here fl 0 = = ,T 0 = ^§?, and 7\ = ^are the 

reflection and transmission coefficients of the two interfaces. The reflection and transmission 

>v 2 2 

coefficients of the layer satisfy the energy conservation relation |3?(o>)| + |T(a>)| = 1 for all co, 
which means that the individual energies of the frequency components of the incoming pulse are 
preserved by the scattering process. The main qualitative difference between the scattering by a single 
interface and the scattering by single layer is that the reflection and transmission coefficients in the 
layer case are frequency-dependent. This frequency dependence originates from interference effects 
between the waves that are scattered back and forth by the two interfaces of the layer. 

2.3. Filtering Property of the Layer 

2.3.1. Reflection 

Let us consider a layer embedded between two homogeneous half-spaces that have the same 
material properties, i.e., the situation in which pi = po and Ki = Kq. We then have R\ = -Rq and 
T\ = To, which implies that the global reflectivity of the layer can be written as: 

l+R%-2Rficos{—) 

The reflectivity is periodic with respect to the angular frequency a> with the period co c = nci/L. As a 
function of the angular frequency the reflectivity goes from the minimal value: 

l^l 2 . = 0 for a) = kco c , k G Z 



'mm 



to the maximal value: 



I^L. = i-lTT-^l for(0 = \k + -J(O c ,ke1 



'max \l + ft^ 

This shows that for any value of the reflection coefficient R 0 of a single interface, there exist 
frequencies that are fully transmitted or fully reflected by the layer. If we consider the case of strong 
scattering T 0 2 « 1, then the transmitted frequency bands have a width of the order of o) c T 0 2 around the 
fully transmitted frequencies kco c . Outside of these bands, where total reflection occurs, the typical 
reflectivity is large, of order 1 — T 0 4 /4. 
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Figure 6. Reflectivity |j?(o>)| versus frequency for a single layer with Rq = -R\ = 0.1 
(a) and R 0 = —R\ = 0.9 (b). The period is co c = nc\IL. 
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2.3.2. Anti-Reflection 

The total transmission phenomenon is also encountered in situations in which the two half-spaces are 
different. Indeed, consideration of human body part as an ideal, fully-transmitting layer is certainly beyond 
perfection. In a microscopic or constituent-wise sense, a human body-part, striated muscle for example, 
iscomposed of water (70.09%), ether-soluble extract (6.60%), crude protein (21.94%), etc. [35]. When one 
is faced with the task of modeling portions of human body as a medium for electromagnetic (signal) 
propagation, there are inevitable practical assumptions and approximations to be made. Thinner layers can 
be considered with more homogeneous characteristics, while for thicker setting with internal variability, the 
aggregate behavior sums up by and large. Thus a certain part of a human body like muscle or fat, when 
considered in macroscopic perspective, can be regarded as a continuum, and hence the idea of a planner- 
layered-medium assumption of human body tissues contextually holds for all practical modeling 
considerations. Such model could greatly affect related system design pertaining to crucial parameters 
(antenna and others); for instance, in addition to the frequency and bandwidth employed for the signal, the 
measurement of the depth of the body at which the transceiver of an implant device would function with 
acceptable accuracy has a lot to do with such parameters as permittivity, permeability, and impedance of 
the intermittent layers of body tissues. 

Figure 7. A propagation system consisting of 3 layers; air, human body channel, and a transceiver. 
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Figure 7 shows a grossly approximated propagation system consisting of three layers: air, human 
body channel, and a transceiver. Admittedly, finding fully-transmitting or fully-reflecting layers in 
body-parts is unrealistic, but it is widely adopted practice to model systems using phantoms comprised 
of components having somewhat homogeneous characteristics, yet closely resembles body-parts in 
aggregate behavior. In some recent works ([33] and [34]) UWB antenna impedance matching has been 
studied in the context of biomedical implants. We assume that the two homogeneous half spaces 
(Figure 7) have different impedances Co ^ C2, then it is possible to choose the thickness L and the 
impedance C,\ of the layer so that a given frequency w will be fully transmitted from one half-space to 
the other one, which would not be the case in absence of such a layer. From the analysis of the 
reflectivity function: 



|#(a))| : 



1 - 



i-rI-rI+rIrI 



(20) 



l+2/?o/?icos(^)+i?gi?f 

One can show that a necessary and sufficient condition for |J2(o>)| to be zero is that 
Ro + R| = — 2R 0 Ri cos (~r~J> I n the case Co ^ Ci this in turn enforces one to choose the impedance of 

the layer to be ^ = ^oh ( so that i?o = Ri) and the thickness L to be chosen so that a>LI(%c\) is half an 
integer (so that cos(2 coLlc\) = -1). Usually the thickness is chosen to be equal to a quarter of the 
wavelength, meaning a>L/(nci) = Vi. 



Figure 8. Transmitivity |T(a>)| versus frequency for a single layer with Rq 
(a) and Rq = ~Ri = 0.9 (b). The period is co c = %C\I2L. 
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2.4. Path Loss in the Human Body (Near Field Far Field Consideration) 



When EM RF waves propagate in freespace, the power received decreases at a rate of (l/d) n , 
n being the coefficient of pathloss. Other kinds of losses would be fading of signals due to multipath 
propagation. However, for propagation of EM waves in a lossy medium like human tissue, the losses 
would be mainly due to absorption of power in the tissue, where it is dissipated as heat. As the tissue 
medium is lossy and mostly consists of water, the EM waves are attenuated considerably before they 
reach the receiver. The Specific Absorption Rate (SAR) is useful in determining the amount of power 
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lost due to heat dissipation. SAR is defined as power absorbed per unit mass of the tissue [29]. SAR is 
a standard measure of how much power is absorbed in the tissue and depends upon E- and //-field 
strengths. By determining the average SAR over the entire mass of the tissue between the transmitter 
and the receiver, we are able to compute the total power lost. SAR in the near field of the transmitting 
antenna depends mainly on the //-field, whereas the SAR in the far field of the transmitting antenna 
depends mainly on the E- field. We use Maxwell's E- and //-fields equations for lossy medium to 
obtain the average SAR of the medium between the transmitting and the receiving antenna in the far 
field and near field, respectively. WBAN applications involve wireless communications between 
implanted biosensor nodes inside human body. 

These nodes exchange data among themselves and also with the base-station. In general, the system 
model consists of numerous biosensor nodes placed inside the various parts of the human body 
surrounded by tissues. In particular, for the development of this model, we consider only one 
transmitting and one receiving antenna separated by a distance d. An elemental short dipole (dipole 
length_wavelength) in a lossy human tissue medium is considered for this purpose [28] and is shown 
in Figure 9. A small area of tissue surrounding the antenna is considered for our analysis. Thus we can 
safely assume the human tissue under consideration to be a homogeneous medium with no sharp 
edges, no rough surfaces and having uniform electric and magnetic properties. The received power is 
assumed to be due only to the power from the transmitter and not from any other source. The space 
around the radiating antenna is divided into near field and far field regions as shown in Figure 10. 

Figure 9. A Hertzian Dipole. 



a 




Figure 10. Field regions around a Hertzian Dipole. 

4 
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The region of space immediately surrounding the antenna is known as the near field region. The 
extent of the near field in the case of short dipoles is given by do = A/2, where X is the wavelength [30]. 
In the near field, the E- and //-field strengths vary rapidly with the distance from the antenna. The far 
field is the entire region beyond the near field. In the far field region, the E- and //-field exhibit a plane 
wave behavior. Power absorbed between the transmitting and receiving antennas can be considered as 
the sum of power absorbed in near field (Pnf) and far field (Pff) regions. The total power absorbed 
between the two antennas is computed by numerical integration. 

Consider an elemental oscillating electric dipole in a lossy medium of conductivity c (S/m), 
permittivity e (F/m), permeability ju (H/m), complex propagation constant y, complex intrinsic 



impedance t] — 



a+ju>e 



[28] at frequency co, as shown in Figure 8. The dipole consists of a short 



conducting wire of length dl, terminated in two small conductive spheres or disks. Assume that the 
current / is uniform and varies sinusoidally with time [28]. The electromagnetic field at a distance 'R' 
for an Hertzian dipole is derived from the vector potential A, given by [28]: 

Hldle- yR 
4n R z z 

where a z is the unit vector in the z-direction, y is the propagation constant, given by y = a + jf> ; 
attenuation constant a and phase constant /? is given by as [28]: 

,1/2 



hie 

a = co — 



J! M=) 2 - 



(Neper/m) 



/» = » = 



i + 



(-)' 



+ 1 



1/2 



(rad/m) 



Spherical components of A (i.e., a R A R + aeA e + a^) are given by A R = A z cos0, A d 
Ap = 0. The magnetic field intensity H and the electric field intensity E is given by [28]: 



-A z sing6 and 



1 1 



E = 



a + jaje 



(V x H) = 



a R 



a + ja)6 L RsinO 39 



Solving the above magnetic and electric field equations for lossy medium and expressing in terms 
of complex impedance rj we get: 



Er=V 

E e =7] 



lldlcosQ 



An 

IdlsinO 
An 



H<p — 



IdlsinO 
An 



R (J_ + M 

\ Y R 3 R 2 ) 

e -yR (JL + L + L) 
Vy/? 3 R 2 RJ 



(21) 
(22) 
(23) 
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2.4.1. Power Absorbed in the Near Field 



The SAR in the near field is given by [31]: 

0 /Id) 

SAR = — — ==(1 + c corr TyHf ms watts /Kg 
P^o 2 + e 2 co 2 

where p is the density of the medium and c cor r is the correction factor to take into account the changed 
reflection properties for small distances R of the antenna from the scatterer. Since we assume both the 
transmitting and receiving antennae are in a same homogeneous medium, the plane wave reflection 
coefficient t is zero. By substituting t = 0 and RMS value of the //-field, the above equation 
reduces to: 

/ . . \ 2 

SAR=-—^==[— 

which gives the value of SAR at a point at distance 'R' and angle from the dipole. Power at 
infinitely small volume (dV — R 2 sin6 dR dO dcp) is: 

AP = SAR x Amass = SAR x p x dV (24) 

The power absorbed in the near field of the lossy tissue can be obtained by computing the average 
SAR over the entire tissue mass in the near field, which is obtained by integrating AP over the entire 
mass in the near field region, i.e., from the surface of the antenna (R = r) to the end of the near-field 
region (R = do): 

P NF = f I I AP 




= a f 



rUQ rTZ r LTZ 

jR=r J 0=0 J(p=C 

ldl\ c r a * r n r 2n . . , _ / 1 \ Y \ 2 2\y\\ 
— R 2 sin 3 6 xe~ 2aR [^ + -^ r + -^-)dRd9d( P 

An J R =r e=o <p=o \R 4 R 2 R 3 



Ver 2 + e 2 a) 2 \ 

Solving by numerical integration and writing -^==^== as |r/|/|y| 



where: 



P NF = <riuof^[A + B + C] (25) 



A = e -2ar (\r\ 2 , d o- r , \Y\Cdo-rT 



2a 4r 2 2r ) 

B c - 2ad J-\Y\ 2 ,dp-r lrl(do-r) \ 
V 2a Ad 2 + 2d 0 ) 

r = e -«(^o + r) ( Xdo-r) 2\ Y \(d 0 -r) \ 
\(d 0 +r) 2 (d 0 +r) J 

The antenna dimensions depend on the wavelength of the wave in the medium given by X m = — [28]. 



Sensors 2010, 10 



5518 



2.4.2. Power Absorbed in the Far Field 

l l 

Neglecting -7,-3 terms from field Equations (21)-(23) for the far field, we have: 



IdlsinO D /V\ 
4n 

IdlsinO D /i^ 

In the far field the specific absorption rate depends only on the E rms value which is given by [29]: 

o 

P 



IdlsinO D /V\ 
IdlsinO D /V\ 



<7 ? 

SAR = —Er ms watts I Kg 



= a -(\ri\\Y\ 



IdlsinO _ aR \ 2 



p\ 4nR ) 

The power absorbed in the infinitely small volume (dV = R 2 sin0 dR dO dcp) in the far field, at a 
distance R and angle 0 from the dipole can again be obtained from (24): 

Idl\ 2 



kP = g(\7]\\y\— ) sin 3 0e- 2aR dRd0 d(p 



The total power absorbed in the far field of the lossy tissue between the source and destination 
antennas can be obtained by computing the average SAR over the entire tissue mass in the far field 
from distance do to d (do is the point where the far field starts). This is obtained by integrating AP over 
the mass in the far field between the two antennas: 

'* r2n I 2 dl 2 dl 



FF 



f f f rdrdl 

AP = a\n\ 2 \ Y \ 2 — ( e -2«do_ e -2« d) 

jR=d 0 Je=oJ<p=o 12na 



Power received 

The effective radiated power (ERP) is obtained by subtracting the loss in the near field (Pnf) and far 
field (P F f between the transmitting and receiving antennas) from the transmitted power P T (i.e., 
(P T -PLoss)G t ), where P Loss = P NF +Pff is obtained from (25) and (26). The power density (P e , Power 
per unit area) at a distance 'd' is different in near field and far field regions: 

(a) P R in the Near Field: There is no general formula for the estimation of field strength in the near 
field zone [30]. Only measurements can provide a simple means of field evaluation. However, 
reasonable calculations can be made for antennas like dipole or monopole. When the receiving 
antenna is in the near field region of the transmitting antenna, the power density does not 
necessarily depend on the distance from the antenna, but varies rapidly with distance, and may 
exhibit oscillatory behavior. The magnitude of on-axis (main beam) power density varies 
according to the location in the near field and its maximum value is approximated by [32] 
P e = 16 8P/nL 2 , where L is the largest dimension of the antenna, P is P T - P NF , and 8 is the 
aperture efficiency (typically 0.5-0.75) [32]. It can be approximated as S = A e /A (A e is the 
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effective aperture and A is the physical area of the antenna). The power received by the 
receiving antenna in the near field can be approximated by: 



Pu = P P A e = 



16S(P T - P NF ) 
nL 2 



A, 



(b) P R in the Far Field: On the other hand when the receiving antenna is in the far field region of 
the transmitting antenna, the power density is dependent on the distance d and is given by: 



Pe = 



(Pi 



Loss 



And' 



G t 



The power received by the receiving antenna in the far field is P R = P e A e , where the receiving 

x 2 

antenna aperture A e is given by A e = — G r Here, G t and G r are the gain of the transmitting and 



G t G r 



receiving antenna, respectively. Thus the received power is: 

p _ (Pi ~ ?nf ~ Pff)^ 2 
R ~ (And) 2 

and a total phase change of is involved during the propagation of the wave. Thus, PMBA can be 
used for calculating the propagation loss using the two Equations (25) and (26). 

Figure 11. PMBA (tissue medium) and Freespace Pathloss at 2.4 GHz. 




3. Numerical Results 



The permittivity of biological tissues depends on the type of tissues (e.g. skin, fat, or muscle), water 
content, temperature, and frequency. However, the permittivity and frequency may also determine how 
far the EM wave penetrates into the body. The term depth of penetration (D p ) usually quantifies this. It 
is observed from Equations (1) and (2) that the wave gets attenuated as it propagates in the biological 
material along the z-axis. As shown in Figure 1, variation of radiation power density has been 
compiled for four different frequencies (27 MHz, 100 MHz, 433 MHz, and 1,500 MHz) with respect to 
the depth in muscle. At a given depth, usage of lower frequency results in a higher power density as 
illustrated in Figure 1. We discovered a distinct feature (demonstrated by Figure 2) which states that 
attenuation is more in fat than that in muscle with respect to frequency. Expressions for left- and 
right- going modes of a pulse have been derived in appendix B. Equations (19) and (20) express the 
reflectivity and transmitivity of a layer respectively. Equations (7) and (8) are the expressions of the 
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left- and right-going modes respectively of a pulse when scattered by a single interface between two 
layers. Figure4. shows the scattering of a Gaussian pulse by an interface separating two homogeneous 
half-spaces (cq, Co, z < 0) and (c\, &, z > 0). The spatial profiles of the velocity field (a) and of the 

>v 2 

pressure field (b) are plotted at times t = -4, t = -3,..., t = 6. Reflectivity | J2(a>) | versus frequency has 
been depicted by Figure 6. We can see here that, at a certain period w c = ncJL, the reflectivity 
maximum which is almost 0.04 for a layer with R 0 = —Ri = 0.1 (a) and almost 1.0 for a layer with 

>v 2 

Ro = ~Ri = 0.9 (b). In Figure 8, Transmitivity |T(a>)| versus frequency curves have been drawn. 
Transmitivity has been found to vary periodically with a certain frequency. The period of this 
frequency band depends upon the choice of the layer thickness L. For a layer with Ro = —R\ = 0.1, 
transmitivity is about 1.0 (a) and for a layer with Ro = -R\ = 0.9, transmitivity is about 0.9 (b). Here the 
period is co c = nc\l2L. From Equation (20), we can see that by the proper choice of the impedance of 
the layer to be ^ = ( so that Ro = Ri) and the thickness L to be chosen so that <x>LI(%c\) is half an 

integer (so that cos(2 coLIci) = -1), we can form a fully transmitting layer. Usually the thickness is 
chosen to be equal to a quarter of the wavelength, meaning coL/(kci) = ¥i. Therefore, from the results 
shown above, we can infer that a layer can either fully reflect or fully transmit any incoming wave at a 
certain frequency or frequency band. We can use these results to UWB by proper choice of impedance 
and the thickness L. Power loss in near filed and far field due to absorption has also been analyzed. A 
propagation loss model (PMBA) for homogeneous tissue bodies has been presented, which compares 

^— ) ; with loss coefficient n = 3). A 

frequency range of (900 MHz to 3 GHz) has been considered here. We have been able to make a 
conclusion that, compared to freespace, there is an additional 30-35 dB of attenuation at small 
distances the far field (Figure 11). This loss increases further with the distance and frequency. It is 
argued that the human body cannot be considered as layered a media. However, we have been able to 
homogenize the human body channel, which is shown in Appendix C. 

4. Conclusions 

Employing UWB in WBAN involves a lot of promise, just as there are a number of relevant 
challenges. We studied the technical feasibility in this regard with a concentration in electromagnetic 
propagation of the signal across human body. Unlike conventional wireless channels, human body 
comes with a great deal of structural complexity requiring significantly different design considerations. 
The reflection and transmission coefficients of human body are heavily dependent upon the dielectric 
constants as well as upon the frequency. In this work, we investigated a layer-wise model for 
electromagnetic propagation across the components of the body in regard to such key aspects as 
scattering, reflectivity, and transmitivity. Naturally, the segmentation in precise layers are not what we 
come across in a body. But, the approximate model employing homogenization could help assess the 
aggregate behavior of the wireless communication involving implant devices, thus guiding the 
potential design issues for antenna characteristics, for instance. We also presented numerical 
depictions of some of the pertinent signal characteristics. From here on, one could expect to further 
improve the model in terms of suitable layering and other parameters of approximations. 
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Appendix A 

Let us assume that a biological medium is infinite in extent, source-free, isotropic, and 
homogeneous. The medium is isotropic if z is a scalar constant, so D (electric displacement field 
vector) and E (electric field vector) are the same in every direction. A homogeneous medium is one for 
which £ (relative permittivity), \x (relative permeability), and a (conductivity) are constants. For this 
case, Maxwell's equations become: 

Vx£ = -— (A.l) 

dt 

VxH = J + (A.2) 

dt 

V-5 = 0 (A.3) 
V • D = 0 (A.4) 

Using following identities: 

VxVx£ = v(v-i?)-V 2 £ (A.5) 
B = juH (A.6) 
D = sE {A.l) 
J = aE (A.8) 

where, e, \i, and c are relative permittivity, relative permeability, and conductivity of the medium, 

H is magnetic field vector, B is magnetic flux density, J is the displacement current. We can find the 
expressions for wave equation: 



r n2 \(-\ 



d d 

V - ua us — 

dt dt 



(A.9) 



In view of the fact that equations governing E and H in the biological material (Maxwell's 
equations) are linear and keeping in mind that any arbitrarily time-varying function can be expressed 
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as a sum of number of sinusoidal functions, time dependence of the fields, E and H, can be given by 
the factor d 0 ' 1 so that: 

d_ 
dt 

d 2 



■ = ja> 



dt 7 



= -co 



Using the relationships in (A.9), the wave equation becomes: 



where: 



V 2 E + y 2 E = Q 



y 2 = co 2 /lis - jcqucT 



= oo jue 0 \ s -j 



COS 



(A. 10) 



(A. 11) 



And where c is the free space velocity (3 x 10 m/s) and y is the propagation constant. This is in 
general, a complex quantity and may be written in the form: 

y = a + j/3 

where the attenuation constant is: 

4ic 



a = 



,1/2 



CO^iS 



1+ 



€ 



+ 1 



(A. 12) 



for — < 1 

s 



a = 



coyjjue 



4i 



K £ J 



(A. 13) 



and the phase constant in radians per meter is: 



4ic 



coss 



1+ 



f " \ 2 

€ 



1/2 



J 



(A. 14) 



for — < 1 

e 



P = COyl JUS 



1 + 0.125 



f " \ 2 
s 



\ £ J 



(A. 15) 



Using Equation (A. 15), the wavelength A, can be determined by: 
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X = 



In 



(A. 16) 



Appendix B 



In this appendix, we use a number of essential transformations of the wave equation that are 
specific to layered media. We consider the particular case in which the parameters of the medium vary 
in a piecewise-constant manner; in other words, we consider a stack of layers made of homogeneous 
media. We study the propagation of a normally incident plane wave, which enables us to reduce the 
problem to the one-dimensional acoustic wave equations [9]. We will see that the problem can be 
recast as a product of matrices corresponding to the scattering of the wave by the successive interfaces 
between the layers. This is a classical setup for waves propagating in this particular type of layered 
media, and it is extremely useful for direct numerical simulations. The equations for the 
three-dimensional velocity u and pressure p are: 



P d i+v P = o 



(B.l) 
(B.2) 



where p is the density of the medium and K the bulk modulus of the medium. These two equations 
correspond respectively to conservation of momentum and mass. The density and bulk modulus are 
assumed to be spatially varying along the z-coordinate. If the initial conditions correspond to a plane 
wave that is normally incident to the layered medium, then the solution of the equations remains 
independent of the transverse variables, and the transverse velocity is zero. The system can then be 
reduced to the one-dimensional wave equations. It should be noted that more general conditions, 
corresponding in particular to point source, require a more general three-dimensional framework. In 
this appendix, we focus our attention to the one-dimensional case. 

In one-dimensional medium the equations for the velocity u and pressure p are: 

(B.3) 



dt dz 
1 dp{t,z) du{t,z) _ 
K(z) dt + dz ~ U 



with p being the density and K the modulus of the medium, which are both functions of the spatial 
coordinate z. We write this system of equations in matrix form: 

d 

~d~z 



pit.z) 




0 


pO) 


d 


p(t, z) 


.u(t, z). 




K(z)- 1 


0 . 


di 


u(t, z). 



A diagonalization of the 2 x 2 matrix gives 

0 p(z) 



o 



= M(z)- 1 



c(z)" 
0 



0 

-dz)- 1 



M(z) 



where: 



M(z) = 



azr 1/2 

-«z)- 1/2 



az) 1/2 



^(z)" 1 ^ 



azy /2 -az) 1/2 
azr 1/2 azr 1/2 
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with c(z) = yjK(z)/ ' p{z) and ^(z) = ^K{z)p{z) being respectively the local speed of sound and 
impedance. The system can then be written as: 



d 

~d~z 



p(t,z) 
u{t, z) 



d 

di 



pit, z) 
u(t, z) 



In this representation the material parameters p and K may either vary or remain constant (special 
case) with respect to the space coordinate z. 



Right- and Left -Going Waves 

We consider the special case with a homogeneous medium in which the coefficients p and K are 
constant. Consequently the speed of sound c and the impedance C, are constant, and the system can 
be written: 



-(M\ P(t ' Z) }) = --\ 1 °]-(M\ Pit ' Z) }) 
dz\ [u{t,z)\) cU> -Vdt\ [u{t,z)\) 



in 



p(t,z)l 



Then if we define: 



A{t,z) 
B(t,z) 



= M 



p{t, z) 
u(t, z) 



C 1/2 p(t,z) + ( 1/2 u(t,z) 
-C 1/2 p(t,z) + t 1/2 u{t,z)_ 



it follows that: 



d 

dz 



A(t,z) 
B{t,z)\ 



= _ i rl 0ia_ 



0 l d \A{t,z) 
B(t,z)\ 



(B.4) 



(B.5) 



The equations for A and B decouple: 

dA{t,z) IdA (t,z) 

dz c dt 
dB{t,z) ld£(t,z) 
c 



= 0 



= 0 



dz c dt 

and the waves can be written A{t,z) = a(t — z/c) and B{t,z) = b(t + z/c) for some wave-shape 
functions a and b. Thus, in the constant medium case we have decomposed the wave into the right- and 
left- going waves A and B, which do not interact. 

To fully specify the problem we have to prescribe initial conditions, for instance the velocity and 
pressure profiles at time t = 0: 

u(t = 0,z) = u 0 (z),p(t = 0,z) = p 0 (z) 
We then translate these initial conditions for u and p into initial conditions for the modes A and B: 

4)(-z) == Ait = 0,z) = r 1/2 PoO) + ( 1/2 u 0 (z) 
B 0 (-z) == B{t = 0,z) = -C 1/2 PoO) + ( 1/2 "oO) 
which gives the expressions for the modes: 

A{t,z) = Aq (ct — z) , 5 (t, z) = S 0 (ct + z) 
and finally the expressions for the wave: 
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, , __ 1/2 i4(t f z)+S(t,z) 
u(t,z) = ( 1/2 

The initial conditions determine the mode decomposition and can be chosen to generate a pure 
right-going wave (if p 0 = £ u o) or a P ure left-going wave (if p 0 = — (u 0 ). 

A more physical way to generate a wave is to assume that the wave vanishes as t -> — oo and to 
introduce a source term in the acoustic wave equations: 

du(t, z) dp(t,z) 



dt dz 



= F(t,z) 



1 dp(t,z) du(t,z) 
K- 1 H \ - + ^ =0 
at oz 

By assuming a point source F(t, z) = Z, 1 ? 2 f(i)8(z), the system for A and B becomes: 

dAjt.z) I dAjt.z) 

+ ~ T— = S(z)f(t) 



dz c dt 
dB(t,z) ldB(t,z) 
dz c dt 



= -5(z)/(t) 



whose solutions are: 



fit - z/c) ifz>0 
0 if z < 0 



0 if z > 0 



4(t,z) = { 

5(t,z) =[f(- t + z/c) ^ z< o 
As a result, the velocity and pressure fields are: 

_ r 1/2 f/(t - z/c) if Z >o 

K ' } 2 {f(t + z/c)ifz<0 

v(t z)=^—i nt - z/c) if2> ° 
VV ' } 2 \-f(t + z/c)ifz< 0 

This means that the source term generates two waves with equal energy that propagates to the right 
and to the left. 



Appendix C 

Homogenization : 



For a channel like the human body where muscle, fat, blood cannot be considered as slabs or layers, 
we consider the idealized situation in which the parameters vary only with depth, and moreover, we 
make the important assumption that the variations are on a relatively fine scale. We assume that the 
scale of variation is small compared to the distance traveled by the pulse, as well as compared to the 
wavelength of the pulse. One may then expect that the waves are not strongly affected by the 
impedance in any particular layer. When a pulse propagates through such fine layers, the interaction 
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with each layer is small, and propagation is not much affected. The pulse therefore travels as if the 
medium were homogeneous with the layers replaced by "averaged" ones. In general, we refer to this 
homogeneous medium as the homogenized medium. It is also referred to as an effective, average, or 
equivalent medium. We start by writing the medium parameters in the formp = p(z/V),K = K(z/V), 
with I a parameter that can be viewed as the layer size. Thus, p(z) is the variable density when 
observed through a magnifying glass with magnification factor l/l. We then observe the fluctuations 
on their natural or intrinsic scale of variation. Typically, we will model p(z) as a stationary random 
process. Again we consider that the slab has constant parameters p = p and K = K, then <f = ^Kp and 
c = sj K/p. Assuming that the local propagation speed c varies with z but the impedance is constant. 
We choose the value ( to be this constant. With these assumptions we have: 

p(z) K c 



We introduce propagator P w for each frequency co, then we have: 

^P w (0,z) = H w (z,z//)P w (0,z), P w (0,0) = / 

where we have denoted the identity matrix by I and we have introduced: 

(AW(z') - 1) A«(z>- 2iwz/c " 



(C.l) 



10) 



-A ( - ) (z')e +2iwz/c " (l - A (+) (z')) 



l/p(z) , K 



A c ±)(z)= _ lT± _ ; 

A< ±> (z)=A( ± '(z/i).z'=j 

The matrix P w (0,z) "propagates" the wave components from z = 0 to any other location z > 0, 
since the linearity of: 



d Ta] _ ia> 
dz [hi c 



(a<+> - 1) 

)g+2iwz/c 



^( - )g-2i<Dz/c 

(l-A«) 



0 < z < L 



(C2) 



implies that: 



a(<x), z) 
b(u), z) 



= ^(0,z) 



a(<x), 0) 
b(a), 0) 



(C.3) 



for any z G [0, L] . Equation (C.2) can be found by Fourier transforming A(t, 0) + P 0 5(t, 0) = 
T 0 f(t) , T 1 A{t,L)-B{t,L) = 0 with boundary conditions a(co, 0) + R 0 b(a), 0) = T 0 f (oj) , 



R x e~ c a(a),L) -b(a),L) = 0, where f(a)) = j f(s)e ia>s ds. 

Equation (C.l) is diagonal and can be integrated by exponentiation: 

r e ik)S;(z) Q 



P w (0,z) = 



0 g-!w5;(z) 



(C4) 
(C.5) 
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To determine the effective medium that emerges in the limit of fine layering Z -> 0, we see from 
(C.5) that we need to study the behavior of J Q Z \ ^~7q ~ ^) dy as Z -> 0. Homogenization can be 

illustrated using the simple model in which the medium is made up of independent and identically 
distributed layers of equal width / -» 0. The medium is defined by one sequence of independent and 
identically distributed positive random variables C n that are bounded and bounded away from zero. 
The local speed of propagation is given by: 

c(z/Z) = C [z/l] 

where [x] denotes the integer part of x. Since [z/l] -» oo as Z -> 0, we can apply the law of large 
numbers to obtain: 

-■= \ dy = I J c~ x {y)dy 




(?) 

= m * M^- 1 t) + ' (f - [*/<])4 - * fc] «*> 

Z 1 . ' 

The convergence is in the almost sure sense, for almost all realizations of the medium, or with 
probability one with respect to the randomness. In this setting, homogenization in the frequency 
domain means that we should choose c such that in the limit that Z -» 0 the propagator P w (0, z) 
becomes the identity for all z. Using (C.4-C.5), we see that we must have: 

( rir" 1 

c 



Thus, the harmonic mean of local propagation speeds is the homogenized or effective propagation 
speed. This effective propagation speed is frequency independent in this example, and therefore it is 
also the effective propagation speed in the time domain. We can get the transmitted and reflected 
waves in the Fourier domain: 

ai(a),L) = 5l/(<o) 

1 + RoRte 21 s 

— Ro + Rie c „ 
b 0 (a>,0)=^ i ^/(a)) 

1 + R 0 R 1 e 21 c 
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